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ABSTRACT 

We investigate through dedicated numerical simulations the evolution of the soft X- 
ray absorption properties of a cloud surrounding a gamma-ray burst source. We show 
that the absorption properties of the material are strongly modified by the ionization 
induced by the intense burst flux. We derive the temporal evolution of the measured 
column density as a function of the density and size of the absorbing medium. Even if 
their statistical significance is not extremely compelling, we find that the detection in 
several bursts of variable absorption during the gamma-ray phase can be accounted for 
if these bursts are associated to overdense regions in molecular clouds with properties 
similar to those of star formation globules. We fit our model variable column density 
to the data of GRB 980329 and GRB 780506, showing that with this method the size, 
density and density distribution of the material surrounding a burst can be measured. 
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1 INTRODUCTION 

An important way to unveil the nature of the Gamma-Ray 
Burst (GRB) progenitors is to investigate the properties of 
the ambient medium surrounding them. The most popular 
models for GRB formation involve either the collapse of a 
single massive star or the coalescence of two compact ob- 
jects such as two neutron stars or a neutron star and a black 
hole (Eichler et al. 1989; Narayan et al. 1992; Woosley 1993) 
These models could be constrained by knowing the GRB en- 
vironment. Massive stars have very short lives, thus they will 
explode in star-forming regions, which are typically charac- 
terized by very dense environments. On the other hand, most 
merging neutron stars would be very old and would have typ- 
ically traveled far from their birthplace. The scenario of com- 
pact merger progenitors could thus be suggested by medium 
to low density environments. Recent observations (see e.g. 
Kulkarni et al. 1999) for the subclass of long bursts seem to 
suggest their association with star-forming regions. 

Within the context of GRB progenitors involving the 
collapse of a massive star, several proposals have been made: 
the so-called "hypernova" scenario (Paczyriski 1998; Mac- 
Fadyen & Woosley 1998) predicts a simultaneous explosion 
of the burst with a Type Ic supernova, while the "Supra- 
nova" model by Vietri & Stella (1998) predicts a two-step 
explosion, in which the burst onset follows a supernova ex- 
plosion by several months. While the rebrightening of the 
optical afterglows detected in several bursts ~ 30 days af- 



ter the burst explosion (Bloom et al. 1999; Reichart 1999; 
Galama et al. 2000) is suggestive of the first scenario, the 
detection of iron emission and absorption features in the X- 
ray emission of the bursts and afterglows (Amati et al. 2000; 
Piro et al. 2000) can be better explained with a two-step ex- 
plosion (Lazzati et al. 2001a). Independent tests of these 
different scenarios are therefore needed. A study of the den- 
sity distribution in the close environment of the bursts can 
be a powerful diagnostic of the type of progenitors. In fact, if 
the burst explodes simultaneously with the supernova, the 
photons propagate through the pre-explosion stellar wind 
[n(r) oc r~ 2 ], while in the two-step scenario a high density 
metal enriched supernova remnant is expected to surround 
the burst explosion site (Lazzati et al. 1999). 

Various methods have been proposed to investigate the 
properties of the GRB environment. It has been shown that 
the effect of the surrounding material on the propagation of 
photons can be a probe of the material itself through the 
detection of absorption features and their time evolution 
(Perna & Loeb 1998; Meszaros & Rees 1998; Bottcher et al. 
1999; Ghisellini et al. 1999; Lazzati et al. 2001abc; Ramirez- 
Ruiz et al. 2001a). In X-rays, it is particularly interesting the 
possibility to detect iron K a absorption edges (Ghisellini et 
al. 1999; Bottcher et al. 1999, Weth et al. 2000) and their 
time evolution (Lazzati et al. 2001c; Bottcher et al. 2001). In 
the optical, the evolution of the equivalent width of absorp- 
tion features has been investigated by Perna & Loeb (1998). 



2 Lazzati & Perna 



On the other hand, modelling the afterglow spectral and 
intensity evolution in the framework of the external shock 
synchrotron model can also constrain the type of environ- 
ment (Wijers & Galama 1999; Berger et al. 2000; Panaitescu 
& Kumar 2001, Ramirez-Ruiz et al. 2001b). In the first case, 
however, the search for the signatures of the material is made 
difficult by the low signal-to-noise ratio of the lines, while 
the second procedure heavily relies on a particular model, 
and hence depends on its assumptions. 

In this paper, rather than considering the feature pro- 
duced by a particular transition, we consider the continuum 
soft X-ray absorption caused by the photoionization of low- 
intermediate atomic number (Z) atoms along the line of 
sight to the burst (Morrison & McCammon 1983). This ab- 
sorption is usually quantified by means of the quantity Nh, 
the column density of neutral material with solar metallicity 
that would cause the observed absorption (in the simulations 
we use abundances from Anders & Ebihara 1982). The pres- 
ence of Nh is generally much easier to detect with respect 
to a Fe K a absorption edge. In fact, the opacity of iron at 
7.1 keV is less than one hundredth of the opacity at 1 keV of 
a material with solar iron abundance (Morrison & McCam- 
mon 1983). On the other hand, the precise measurement of 
the absolute Nh value can be model dependent. In fact, the 
Nh measurement is performed by computing the difference 
of the observed soft X-ray flux with respect to an assumed 
model in the same band. A spectral break at ~ 1 — 2 keV may 
for this reason induce a high spurious Nh measurement. 

In the following sections we present the results of nu- 
merical simulations for the evolution of the measured Nh 
as a function of the time elapsed after the onset of a source 
of ionizing photons in its centre. We consider in particular 
the possibility of detecting a variation of the column density 
during the ~ 100 seconds of duration of a typical GRB. By 
imposing that the Nh should remain above the detection 
threshold for a reasonable time (Is) but vary in less that 
the GRB duration, we can strongly constrain the geomet- 
rical properties of the material surrounding the GRB. We 
then compare our results with the Nh measurements per- 
formed in the early emission of GRB 980329 (Frontera et al. 
2000) and GRB 780506 (Connors & Hueter 1998). 



2 NUMERICAL SIMULATIONS 

In the case of a single atomic transition it is possible to de- 
rive analytical equations to describe the time evolution of its 
opacity in the optically thin regime (Lazzati et al. 2001c). 
In the case of continuum absorption the calculation is com- 
plicated by the need of keeping track of all the ionization 
states of the various elements and relative cross sections. 
For this reason the problem can be faced only with the use 
of numerical simulations. 

We have used the time-dependent photoionization code 
described in Lazzati et al. (2001c), which is based on routines 
developed by Raymond (1979) and specialized to the case 
of GRBs and their remnants by Perna, Raymond & Loeb 
(2000) and Perna & Raymond (2000). As described by Perna 
et al. (2000), this code computes in a fully time-dependent 
fashion the ionization state of the 12 most abundant astro- 
physical elements, while updating, at each time step, the 
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Figure 1. Opacity in the range [2-10] keV for a cloud with solar 
metallicity, _R = 3 X 10 18 cm and initial column density Nh(0) = 
3 X 10 21 cm -2 . In the main panel, from top to bottom, we plot 
the absorption at times t = 0, 10, 40, 80, 120 and 160 seconds. 
In the inset, the column density is shown as a function of time. 
Filled dots underline the column densities for which a spectrum 
is plotted in the main panel. 

temperature of the plasma^]. The ionizing continuum is sim- 
ulated in the range [0.01-10] keV with a spectrum relevant 
to GRB emissionn. We adopted a photon index F = 1 and 
constant luminosity Lx = 1-3 x 10 49 erg s _1 in the [2-10] 
keV band for t > to. We considered three different radial 
distributions of the density of absorbers: a uniform cloud 
with R = 3 x 10 18 cm, a shell with R = 3 x 10 18 cm and 
AR/R = 0.1 and a wind with n(r) = no(R/ Rmin) 2 for 
R > Rmm = 10 13 cm and n(r) = for smaller radii. For 
the uniform and shell geometries, we performed simulations 
with initial column densities ranging from 3 x 10 19 cm -2 
to 3 x 10 24 cm -2 , logarithmically spaced by a factor vTo. 
In the case of the wind geometry, we considered only high 
values of the initial column density because of its very rapid 
decrease with time, and we performed simulations for the 
initial values Nh (0) = 10 22 , 10 23 and 10 24 cm -2 . 

Even though the simulations for the uniform and shell 
geometries are performed for a single value of the outer ra- 
dius of the absorber distribution (see above), these can be 
applied to rescaled matter distributions. In fact, given a dis- 
tribution of absorbers n(r) and the ionizing luminosity Lx, 
the results of the simulations can be extended to any dis- 
tribution n'(r) = a" 1 n(r/a) (where a is a scale parameter) 
and any luminosity L' x ~ bLx- The time of photoioniza- 
tion of an ion is proportional to the flux, which scales as 
F oc L/r 2 , where r is the distance of the ion from the pho- 

1 See also Schwarz 1973 for related work. 

2 Note that the ionizing continuum has to be simulated in a spec- 
tral band wider than the one in which the absorption is measured 
since soft photons can ionize ions that contribute to the opacity 
at higher frequencies. 
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Figure 2. Evolution of the column density with time for a uni- 
form (solid line) and shell (dashed line) environments. In both 



cases the initial column density is Nu(0) 
outer radius is R = 1 pc. 
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ton source and L the luminosity of the ionizing photons. For 
this reason 



N H{aM {t) = N H (^) 



(1) 



where Nn(t) is the column density as a function of time 
from the simulation and Nn^ a b y (t) is the column density as 
a function of time for the ionizing luminosity L' x = bLx 
and the distribution of absorbers n'(r) = a -1 n(r/a). 



software for X-ray data analysis (XSPEC; Arnaud 1996) ne- 
glects the change of cr(v) with ionization to compute the 
column density (model WABS) when the single ion features 
cannot be resolved^. The consequence of this assumption is 
that the column density at a given time has a small depen- 
dence on the band in which it is measured. In the present 
paper we consider the energy range [2-10] keV, since our aim 
is to compare our results with the measurements performed 
with the BeppoSAX WFC. 

In the above approximation, we compute the column 
density as: 



N H (t)=N H (0) 



r(y,t) \ 
r M)/ [2 - 10 ] 



(2) 



where the symbol Op-io] represents the average over the 
frequency range [2-10] keV. In Fig. ^ we show the results of 
a simulation for a uniform cloud with Nh(0) = 3x 10 21 cm -2 
and R — 3x 10 18 cm. The opacities as a function of frequency 
are plotted for several times after the burst onset, while the 
inset shows the evolution of the column density in the first 
100 s of the simulation. 

In Fig. y we show the results of two simulations with 
identical initial conditions, but for different geometries of the 
absorbing material. The evolution of the absorbing column 
density in a uniform cloud and a shell are shown with a 
solid and a dashed line, respectively. The column density of 
a shell environment has a roughly linear decay, slower than 
that of the uniform cloud. This is due to the fact that in a 
shell geometry all the material is located at a large distance 
from the source of ionizing photons, and hence the ionization 
time is longer. With a good dataset, this difference can be 
succesfully used to investigate the radial structure of the 
absorbing medium. 



2.1 Computation of the "observed" column 
density evolution 

The above described numerical code produces frequency and 
time-dependent opacities cr(v,t) by fully accounting for the 
spectral variation with time of the opacity as the gas makes 
the transition from neutral to fully ionized; as such, however, 
the output of the code contains much more information than 
the observational data. This is particularly true for the GRB 
prompt X-ray emission, usually detected with proportional 
counter detectors, such as the BeppoSAX Wide Field Cam- 
eras (WFC). In particular, what is usually measured from 
the observed spectra is the amount of soft X-ray absorption, 
parameterized with the hydrogen column density in the as- 
sumption that all the material is cold, i.e. neutral (Morrison 
& McCammon 1983). We hence applied an analogous pro- 
cedure to our synthetic dataset, averaging the opacity in 
the observed band in order to reduce the information con- 
tained in the simulated spectrum to a column density to be 
compared with the observations. This is equivalent to as- 
sume that the time and frequency-dependent opacity can 
be decomposed as r(u,i) = Nn(t) a(v), where <j(v) is the 
average cross section at frequency u, independent of time, 
weighted with the abundance of each element. Our proce- 
dure is also justified by the fact that the most widely used 



3 BURST ENVIRONMENT 

In several bursts, detection of variable low-energy X-ray ab- 
sorption has been claimed (GRB 780506, Connors & Hueter 
1998; GRB 980329, Frontera et al. 2000; GRB 990705, Am- 
ati et al. 2000; GRB 010222, in't Zand et al. 2001). In this 
section we examine under which conditions such an evolu- 
tion of the measured column density can be detected during 
the prompt emission of the bursts. 

In order to detect an Nh variation it is sufficient (but 
also necessary) that the following two constraints be sat- 
isfied. First, the spectrum of the GRB must be integrated 
in order to perform a measurement of column density. The 
minimum integration time (see e.g. Frontera et al. 2000) is 
of the order of several seconds. Variations of the column 
density can then be detected only if 



3 Note that there are cases, such as that of the warm absorbers, 
where the effects of the spectral variations of the opacity can 
actually be observed through the enhancement of the absorption 
edges of intermediate Z elements (see, e.g., Fabian et al. 1994 for a 
more refined treatment of warm absorbers). In our case, however, 
the situation is made more difficult by the fact that the absorber 
is not in the same ionization state at a given time t. The inner 
regions will be completely ionized at a time in which the outer 
regions are still cold, while gas at intermediate distance is warm. 
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Figure 3. Contour plot of Nh{1) (solid contours) and JVjf (100) 
(dashed contours) in the radius vs. column density parameter 
space. From lower left to upper right, the solid contours show 
JVjj-(l) = 10 21 , 10 22 and 10 23 cm -2 respectively. The same lev- 
els are plotted as dashed lines for 7V//(100). The shaded region 
shows the area in which both conditions expressed by Eqs. (^) 
and (jij) are satisfied. The line-shaded area shows the region in 
which the [2-10] keV band is obscured for the whole burst dura- 
tion (JVh(100) > 10 23 cm" 2 ). 



N H (1) > N Hg ^ ~ 10 2 



(3) 



where Nh(1) is the measured Nh one second after the burst 
onset, and Nh GB/1 is the Galactic column density. With this 
condition we ensure that the column density in the surround- 
ings of the GRB gives a detectable signal in the first inte- 
gration bin. 

On the other hand, in order to detect a variation, the 
column density must sensibly decrease during the GRB 
emission time. Our second condition is then that the ob- 
served column density decrease by a factor of two in less 
than 100 seconds, the average burst duration! 4 : 



N H (100) < N H (l)/2. 



(4) 



We must also consider that the burst itself has to be 
detected. Since the soft 7-rays, in which GRBs are usually 
detected, are not attenuated by phot oabsorpt ion, we con- 
sider Thomson scattering as the most effective attenuation 
process. We then have N H (Q) < l/cr T ~ 1.5 x 10 24 cm" 2 . If 
N H (0) £ 10 23 cm" 2 , the [2-10] keV flux is completely ab- 
sorbed and only a lower limit to the column density can be 
measured. In the extreme case of a large and dense cloud 
(line-shaded region in Fig ^, |i| and []) the opacity in the [2- 



4 The factor 2 that we use in Eq. |ij has been chosen since a siz- 
able decrease of the absorbing column is needed in order to have 
a significant measurement of its variation (see below). We tested 
different values (namely 1.1 and 10) finding that the allowed re- 
gion changes slightly, but the conclusions are not significantly 
affected. 
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Figure 4. Same as Fig. Q but for a shell geometry. The shell 
width is assumed to be 10 per cent of its radius. 



10] keV range is larger than unity for the whole duration of 
the bursts, and no flux can be detected in this band. 

In Fig. ^, ^ and |E| we show that conditions (^) and Q 
can be satisfied only in a limited range of cloud radii, ge- 
ometries and initial column densities. In Fig. [] we show the 
column density at time t — 1 s and t = 100 s as a function of 
the cloud radius and initial column density with solid and 
dashed contours, respectively. The shaded area underlines 
the region in which both conditions are satisfied. This shows 
that, in order to detect Nh variations during the prompt 
emission of the burst, the absorbing region surrounding the 
GRB must be compact, with a radius R < 5 pc and initial 
column density Nh(0) > 10 21 cm" 2 . This translates to vol- 
ume densities n ^ 10 3 cm" 3 (dotted lines). Figure ^ shows 
the same results but for a shell with AR/R — 0.1. 

The possible contribution of a wind environment to the 
column density evolution is considered in Fig. [| In the fig- 
ure we plot the measured Nh as a function of time for 
three different initial column densities: from a thin Nh(0) = 
10 22 cm" 2 wind to a maximum value of Nh(0) = 10 24 cm" 2 . 
Figure ^ shows how, no matter the initial value of the col- 
umn density, it drops to an undetectable level in less then a 
fraction of a second. This is due to the fact that, in the case 
of a wind geometry, the column density is dominated by the 
region very close to -Rmin, for which the ionization time is 
very short. 

In order to compare the inferred cloud properties with 
the properties of Galactic molecular clouds, in Fig. ^ we 
compare the gray region in the parameter space with the 
properties of Galactic molecular clouds and their denser re- 
gions. We assume that the burst explodes in the centre of 
these regions. In Fig. ^, the circles show the properties of a 
sample of cloud radii and masses presented in Leisawitz et 
al. (1989), from which we derived average column densities. 
There are no circles that fall in the allowed region shaded in 
grey. Inside molecular clouds, however, there are overdense 
regions, where star formation is supposed to take place. The 
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Figure 5. Measured column density as a function of time in 
a wind geometry with i? m in = 10 13 cm. The three lines show 
different initial column densities (see inset). In all cases, after a 
fraction of a second the observable column density drops to values 
well below the detection limit. 



triangles show the position in the diagram of a sample of 
dense molecular cloud cores, selected for their presence of 
H2O masers and in which massive star formation is taking 
place (Plume et al. 1997). Their radii and column densi- 
ties are consistent with the shaded region. In some cases 
the column density is even too high. A sample of Bok glob- 
ules (Launhardt et al. 1998) , thought to be directely related 
to star formation, is shown with asterisks. Their properties 
are fully consistent with the shaded region. Finally, we plot 
with diamonds the properties of dense cores in the Taurus 
molecular clouds, where low-mass star formation is though 
to happen (Zhou et al. 1994). These are too compact with 
respect to the shaded regions. It is however worth mention- 
ing that the properties of these overdense regions at the time 
of the burst explosion may be somewhat different, due to the 
radiation pressure exerted by the luminosity of the massive 
star progenitor to the burst. 

Should these Nh variations be confirmed, they indicate 
that the burst progenitors (or at least some of them) are 
located inside overdense regions of molecular clouds. The 
properties of these regions are consistent with those of the 
places in which massive star formation takes place. 



3.1 GRB 980329 

To date, three GRBs detected by the WFC on board 
BeppoSAX showed evidence of variable soft X-ray absorp- 
tion. Published data are available for GRB 980329 (Frontera 
et al. 2000). These data are shown, along with their uncer- 
tainties, in Fig. [?]. The error bars on the last two measure- 
ments, non quoted in the original paper, have been assumed 
equal to the measurements themselves. The timescale has 
been reported to the beginning of the flaring activity of the 
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Figure 6. Comparison between the radius and column densi- 
ties inferred from Fig. jj| and the properties of Galactic molecular 
clouds and some hierarchical structures embedded in them. The 
gray and line shaded areas have the same meaning of Fig. ^, i.e. 
the region in which a column density variation should be observed 
during the prompt GRB phase. Circles show radii and average 
column densities of a sample of Galactic molecular clouds (Lei- 
sawitz et al. 1989). Asterisks refer to Bok globules (Launhardt et 
al. 1998), diamonds refer to dense cores in the Taurus molecu- 
lar cloud (Zhou et al. 1994) and triangles refer to massive cloud 
cores (Plume et al. 1997). The two filled dots are the best fit val- 
ues to the N p measurements of GRB 780506 and GRB 980329 
(see Sect. 3.1 and 3.2). 



burst, which starts ~ 23 s after the trigger (see Fig. 1 of 
Frontera et al. 2000) 

We fitted our uniform-cloud and shell models for the 
evolution of the column density to the data. For the uniform 
model we obtained the best values Nh(0) = 1.8 x 10 23 cm~ 2 
and r — 0.13 pc (n ~ 4.5 x 10 5 cm -3 ), while in the shell case 
we obtain N H (0) = 1.8 x 10 23 cm" 2 and r = 0.066 pc. In 
Fig. ^| the best fit models are shown with a solid line (uni- 
form) and dashed line (shell). The shell model has a better 
\ 2 value, but due to the large uncertainties in the measure- 
ment, the two geometries cannot be significantly disentan- 
gled. The upper right inset shows the la, 90 per cent and 
99 per cent confidence contours for the two parameters in 
both cases. Since the statistical significance of the variation 
of the column density is not compelling, the 99 per cent con- 
fidence contour is not closed. The different evolution of the 
two models shows, however, that with good quality data it 
is in principle possible to understand what is the geometry 
of the medium surrounding GRBs. 

A problem in the fit of GRB 980329 data is the lack 
of a redshift measurement of the burst. In fact all our 
simulations are performed in the rest frame. The model 
should be however modified in two ways if the source of 
absorption lies at redshift z > 0. First, the timescale is 
stretched by a factor (1 + z). Most importantly, if a col- 
umn density Nh„ surrounds a source at redshift z, the soft 
X-ray absorption measurement will yield a measured value 
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Figure 7. Evolution of the column density measured in the 
prompt X-ray spectrum of GRB 980329 (Frontera et al. 2000). A 
uniform and shell geometry for the absorbing column have been 
tested. The best fit models, integrated over the same time inter- 
vals of teh data, are shown with a solid and a dashed line, for the 
uniform and shell geometries, respectively. The inset shows the 
lcr, 90 per cent and 99 per cent confidence contours for the two 
fitted parameters. Again, solid contours refer to uniform density 
and dashed contuors to a shell geometry. 



Nh = (1 + z)~ 2 ' 5 Nh z , where the exponent —2.5 holds in 
the redshift range < z < 4. This implies that for a source 
at redshift z — 1 the measured column density is underes- 
timated by a factor ~ 5. The modelling of the data with a 
known redshift may hence give a best fit which is an even 
more compact and dense cloud. 



3.2 GRB 780506 

Besides the three bursts detected by the BeppoSAK WFC, 
a measurement of column density was possible in a burst 
serendipitously detected in 1987 by the HEAO-1 satellite 
with the A-4 and A- 2 instruments (GRB 780506, Connors 
& Hueter 1998). The burst spectrum was fitted with a three 
component model: a power-law, a thermal component and 
an absorption component, with the column density free to 
vary. The lightcurve was integrated in time bins of 10.24 s 
each and the six spectra were fitted indipendently. The mea- 
sured Nh values are shown in Fig. |j| In the first time bin 
a very significant (6<r) soft X-ray absorption was detected, 
while the fits in the remaining 5 time intervals yielded only 
upper limits. Fitting our uniform cloud model to the data, 
again without the possibility to correct for the redshift of 
the burster, we find a best fit cloud with R = 2 x 10 17 cm 
and Nh{0) = 2 x 10 23 cm -2 . Since only one positive mea- 
surement was possible, as shown by the inset in Fig. ^| the 
radius and initial column density of the cloud cannot be pre- 
cisely determined simultaneously. We can safely conclude, 
however, that the density of the cloud was n Si 10 6 cm -3 . 
Since for this burst only one positive measurement was ob- 



22 



10 



21 




N H (0) (enr 2 ) 



10 



20 30 
time (s) 



10 



50 



Figure 8. Same as Fig. | but for GRB 780506 (Connors & Hueter 
1998). The solid line shows the best fit model with N H (0) = 
2 X 10 23 cm -2 and R = 2 X 10 17 cm integrated over the observed 
time intervals. The dashed line shows the best fit model without 
the time integration. The inset shows the lcr, 90 per cent and 99 
per cent contour levels for the two parameters and, with dashed 
lines, isodensity contours for n = 10 6 and 10 7 cm -3 . 



tained, it is not possible to distinguish between a uniform 
and a shell geometry. 



4 DISCUSSION AND CONCLUSIONS 

We have performed numerical simulations to compute the 
evolution of the soft X-ray absorption of a cloud under the 
ionization flux of a Gamma-Ray Burst exploded at its centre. 
We have considered a uniform cloud, a shell and a wind 
geometry, and we have shown that variations of the column 
density during the duration of a GRB can be detected only 
under well-defined geometrical conditions of the absorbing 
material. 

The column density of a wind, in which the bulk of the 
material responsible for the absorption is close to the GRB, 
cannot survive the first milliseconds of the illumination and 
is undetectable with the present instrumentation. On the 
other hand, if the material is spread uniformly over a region 
with size 0.05 R 5 pc, the soft X-ray absorption can 
be observed to evolve in time. Present data do not allow 
to firmly detect the evolution of the column density, nor to 
tell if the material is uniformly spread around the source or 
rather concentrated in a shell at distance R from the burster. 
By comparing the allowed properties of molecular clouds 
with the results presented here, we conclude that the weak 
evidence of Nh variation during the prompt X-ray emission 
of some GRBs can be accounted for if they are associated 
to compact overdense regions similar to the cocoon of star 
formation within molecular clouds or to Bok globules. 

A number of approximations have been applied to 
this analysis and are worth a discussion. First, we con- 
sider an ionizing flux with fixed spectrum and luminos- 
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ity. Both these assumptions are coarse approximations of a 
GRB lightcurve and spectral evolution. Even though GRB 
are frantically variable, the results of our simulations are 
exact if the time t is substituted with the quantity if — 
f Q L(t) At/ J ' GRB L(t) dt, where t G RB is the total burst du- 
ration. Since in this work the detailed time variability of 
the Nh lightcurve is not addressed, the approximation of 
constant luminosity is reasonable. 

A more important issue is the variation of the spec- 
tral slope with time. More complete simulations of time- 
dependent column densities will be presented in a forth- 
coming paper (Perna et al. in preparation), in which the 
dependence of the observed column density on the ionizing 
spectrum will be discussed in detail. However, F = 1 is a 
good value for the average time integrated low-energy slope 
of GRBs (Preece et al. 2000). 

We have assumed that the medium is at rest during the 
photoionization process. While Lazzati et al. (2001c) showed 
that the fireball interaction is not important, radiation pres- 
sure may accelerate the medium. The ratio of the timescale 
of photoionization t p i over the timescale of acceleration t acc 
is given, for aT = l spectrum, by: 

tpi _ hv h a T 6 / hut W hvth \ 2 3 

tacc m p c 2 a pi ~ UllkeW VlkeV J { ' 

where v\, is the break energy of the GRB spectrum and 
o-pi ~ 2 x 1Q~ 22 (huth/i keV) -2 ' 3 is the average photoion- 
ization cross section at the threshold energy /li^h weighted 
with the element solar abundance (Anders & Ebihara 1982) . 
Even for iron, the relevant element with the highest thresh- 
old energy, the acceleration timescale is several orders of 
magnitude longer than the photoionization one. 

Finally, as discussed in §4.1, the redshift of the burst 
introduces time dilation and decreases the observed column 
density at a given time. All these factors have not been taken 
into account in Fig. ^ and Fig. ^, so that the gray shaded 
regions should be considered indicative, with uncertainties 
up to a factor of order unity. Nevertheless we can firmly 
conclude that if a column density variation is observed in the 
first ~ 100 seconds of the prompt emission of a GRB, then 
it must be located within a cloud with properties roughly 
consistent with those of Fig. || or inside a shell (Fig. |^). 
We have not considered here the possibility of a beamed 
fireball. In the case of a conical jet, as discussed to date for 
GRB fireballs (Sari et al. 1999, Ghisellini & Lazzati 1999, 
Panaitescu & Kumar 2001), beaming does not affect our 
calculations, since all processes take place along the line of 
sight. In the case of a cylindrical jet, the radial dilution of 
photons would be different, and our computations are not 
appplicable. 

Variations in the soft X-ray absorption have indeed 
been claimed in several bursts. Only one of them, however, 
shows this evidence quite clearly. GRB 780506, serendipi- 
tously observed by the HEAO-1 satellite (Connors & Hueter 
1998) has a highly significant soft X-ray absorption in the 
first time integration bin, which vanishes in the rest of the 
burst lightcurve (Fig. g). GRB 980329 (Frontera et al. 2000), 
shows a very well defined trend of the Nh to decrease with 
time, which can be fit in our model (see Fig. Q), but with 
large statistical uncertainties. GRB 990705 (Amati et al. 
2000) shows evidence of time dependent iron absorption fea- 
tures. An evolution of the column density is also claimed, 



but the data are not published. Finally, GRB 010222 (in 't 
Zand et al. 2001) shows clear evidence of Nh variations, but 
the column density appears to rise and fall, a behavior that 
cannot be explained in our model and is likely due to the 
presence of a spectral break that moves inside and outside 
the observing window. 

As a consequence of the high density (n > 10 3 ) de- 
scribed above, the afterglow of the bursts should evolve more 
rapidly than in a very low density environment. In partic- 
ular, the time at which the relativistic to non-relativistic 
transition takes place is t nr ~ 10 {E^/na) 1 ^ d (see, e.g. 
Wijers, Rees & Meszaros 1997), where £53 is the isotropic 
fireball energy in units of 10 53 erg and ne the density of the 
surrounding medium in units of 10 6 cm -3 . Unfortunately, 
there are not accurate afterglow lightcurves for any of the 
bursts described in this paper. Indeed, a very high density 
of the surrounding medium is required in GRB 990705 (Laz- 
zati et al. 2001a). In the case of GRB 010222, which displays 
Nh variations but cannot be straightforwardly included in 
our model, a very high density of the surrounding medium 
is claimed (in 't Zand et al. 2001) on the base of spectral 
and lightcurve considerations. 

The HETE2 satellite has on-board two experiment, the 
Soft X-ray Camera (SXC) and the Wide Field X-ray Moni- 
tor (WXM) potentially sensitive to Nh variations. The SXC, 
has an effective area of ~ 10 cm 2 , is sensitive in the [0.5- 
10] keV range and has a field of view of 0.91 sr. The WXM 
has a larger effective area (175 cm 2 ) and field of view (1.6 sr) 
but is sensitive in the harder [2-25] keV band, where con- 
tinuum absorption is less severe. The combination of these 
two instruments, however, will hopefully clarify the proper- 
ties and time variability of the soft X-ray absorption in the 
early prompt X-ray spectra of GRBs. 
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